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Abstract: Headrick and Takayanagi showed that the Ryu-Takayanagi holographic entan- 
glement entropy formula generally obeys the strong subadditivity (SSA) inequality, a fun- 
damental property of entropy. However, the Ryu-Takayanagi formula only applies when the 
bulk spacetime is static. It is not known whether the covariant generalization proposed by 
Hubeny, Rangamani, and Takayanagi (HRT) also obeys SSA. We investigate this question 
in three-dimensional AdS-Vaidya spacetimes, finding that SSA is obeyed as long as the bulk 
spacetime satisfies the null energy condition. This provides strong support for the validity of 
the HRT formula. 
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1. Introduction 

Strong subadditivity (SSA) is a fundamental property of entropy and a cornerstone of quan- 
tum information theory. It states that, given a tripartite quantum system ABC and a joint 
density matrix p{ABC), the entropies of the subsystems obey the following inequalities: 

S{AB) + S{BC) - S{ABC) - S{B) > (1.1) 
S{AB) + S{BC)-S{A)-SiC) >0. (1.2) 

SSA is a general theorem that depends only on basic facts about Hilbert spaces and the 
definition of the von Neumann entropy. Several proofs of it exist but none is very simple. 
An elementary consequence of SSA is that the mutual information I{A : B) := S{A) + 
S{B) — S{AB), which measures the total amount of correlation between the subsystems, 
increases when we adjoin another subsystem to either argument: I{A : BC) > I{A : B). As 
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a monotonicity property — in space, rather than time — SSA is thus in some sense analogous 
to the second law of thermodynamics, and it plays a similarly fundamental role in quantum 
information theory to that played by the second law in statistical mechanics. 

It was shown by Headrick and Takayanagi [1] that Ryu and Takayanagi's (RT) conjec- 
tured holographic entanglement entropy formula [2, 3] generally obeys the SSA property. (See 
Section 2 for a brief review of the RT formula and this proof.) As no proof or derivation of 
the RT formula exists, the fact that it obeys this highly non-trivial property constitutes a 
crucial check on its validity. The proof is remarkably simple — far simpler than any of the 
proofs of SSA itself. The fact that the RT formula obeys SSA is analogous to the black-hole 
area theorem (the event horizon area increases with time) , which was a crucial motivation for 
Bekenstein and Hawking's identification of the horizon area with the black hole's thermody- 
namic entropy. 

The RT formula is restricted to states represented holographically by static spacetimes 
and to regions of the boundary lying inside constant-time surfaces. Hubeny, Rangamani, 
and Takayanagi (HRT) [4] subsequently proposed a covariant generalization, in which the 
entanglement entropy is given by the area of a certain extremal codimension-2 spacelike 
surface (reviewed in Section 2). While the HRT formula has been applied to a variety of 
time-dependent systems, it has so far been subjected to somewhat less stringent testing than 
the original RT formula. In particular, Headrick and Takayanagi's proof of SSA for RT cannot 
be straightforwardly generalized to HRT, and it is far from obvious whether or why the areas 
of extremal spacelike surfaces obey these global inequalities. Deciding whether HRT obeys 
SSA is thus an important question, for two reasons: First, it is a strong test of that very useful 
formula. Second, if true, the statement that HRT obeys SSA would be a very interesting new 
theorem in general relativity, which — much like the black-hole area theorem — would give new 
insight into the relationship between areas and entropies in quantum gravity. 

In this paper we begin to address the question of whether HRT obeys SSA by a study of 
examples. We restrict ourselves to the simplest non-trivial context for field-theory entangle- 
ment entropies, namely single intervals in two-dimensional field theories. (In order that not 
only A, B, C be single intervals but also their unions AB, BC, ABC appearing in (1.1), (1.2), 
we take them to be adjacent and in the order A-B-C.) According to the HRT formula, the 
entropy of such an interval is simply given by the length of the shortest spacelike geodesic in 
the three-dimensional bulk connecting its endpoints (subject to a certain homology condition 
[5]).^ Despite its relative simplicity, this context allows us to address the essential issues 
that arise in going from the RT to the HRT formula, including having a non-static bulk and 
having boundary regions that do not lie on a constant-time surface. More specifically, we 
consider planar thin-shell AdSa-Vaidya spacetimes, as these are the simplest time-dependent 
asymptotically AdSs exact solutions to the Einstein equation. In particular, their simplicity 
allows us to solve the geodesic equation exactly (using the method of [6]), and systematically 

'^Because of the homology condition, in certain cases the extremal surface may also include other connected 
components that do not extend to the boundary (for example that wrap the horizon of a black hole). In the 
present investigation, for simplicity we will neglect this possibility. 
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explore the space of spacelike geodesies. ^'^ In order to investigate the role of the null energy 
condition (NEC), we consider both positive- and negative-energy infalling matter. In Section 
3, we study intervals that lie on a constant-time surface of the boundary (with respect to the 
Lorentz frame distinguished by the system's spatial translation symmetry). In Section 4, we 
turn to more general spacelike intervals. This is an important case to consider when hunting 
for violations of SSA, because, as we explain there, certain such configurations saturate SSA 
in the static limit. 

Our investigation of Vaidya found no violations for positive-energy infalling matter. If 
violations of SSA were generic, we would expect them to appear in the Vaidya case, so 
these results provide a highly non-trivial check on the HRT formula. On the other hand, we 
observed that negative-energy matter leads to violations, which indicates that, if the HRT 
formula is correct, the NEC somehow plays a crucial role in maintaining its consistency. (The 
correlation between SSA and the NEC was also noted by Allais and Tonni [10], based on the 
study of constant-time intervals.) The NEC plays a key role in assuring the consistency of 
the area-entropy relationship in other gravitational contexts, including the area theorem and 
the covariant holographic entropy bound [13, 14], so its appearance here is very natural.^ 

While our positive results for AdSa- Vaidya are a good start, there are clearly many more 
cases that should be checked in order to be confident that HRT obeys SSA generally. An 
obvious generalization would be to higher dimensions. Another would be to more complicated 
shapes and configurations of boundary regions. A third would be the case where the bulk has 
non-trivial topology, and the homology condition comes into play. Fourth, there are many 
cases (with and without non-trivial bulk topology) where there is a phase transition and the 
minimal extremal surface switches from one topology to another as the boundary region is 
varied. Finally, in order to better understand the role of the NEC, one could investigate in 
examples whether there is a direct connection between the violation of SSA and the violation 
of holographic c-theorems that are seen when the NEC is violated [15, 16]. 

Of course, ultimately one does not want merely a laundry list of examples where SSA 
is satisfied (assuming it is indeed satisfied). Rather, one wants to understand how the HRT 
formula pulls off such a trick: one wants a proof, perhaps along the lines of the proof of the 
area theorem. As argued above, such a proof would not only provide very strong support to 
the HRT formula, but would almost certainly provide new, deep insight into the connection 
between areas and entropies in gravitational theories, a connection about which we have so 
many hints yet still so little real understanding. 



■^The same investigation was also carried out for thick-shell Vaidya spacetimes, yielding substantially identi- 
cal results. Here we report the thin-shell results, since those were obtained analytically whereas the thick-shell 
case required numerical integration of the geodesic equation. 

^Several papers, including [7, 8, 9, 6, 10, 11, 12], have previously considered entanglement entropies in 
AdSa- Vaidya spacetimes (mainly focusing on the case of constant-time intervals). 

''Recent work on holographic c-theorems based on entanglement [15, 16] also indicated that the NEC is a 
necessary condition. 
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2. Review: Holographic entanglement entropy formulas 



2.1 RT 

Let us briefly recall the Ryu-Takayanagi (RT) formula [2, 3], and the proof that it obeys SSA 
[1]. The formula applies to field theories that are holographically dual to Einstein gravity, 
in states that are described by static classical geometries. Staticity implies the existence of 
distinguished constant-time slices, both in the bulk and on the boundary; the RT formula 
gives the entanglement entropy of subregions of such a boundary slice, and we work entirely 
within the corresponding bulk slice, with its Euclidean induced metric — the time direction 
and components of the metric play no role. Specifically, given a boundary region A, the RT 
formula says that its entanglement entropy is given by the area of the minimal surface in the 
bulk within a certain topological class:^ 

S{A) = min area(m) ; (2.1) 

m ^ A means that m and A are homologous, i.e. there exists a bulk region r such that 
dr = myj A (and therefore necessarily dm = dA, i.e. m is "anchored" on the boundary along 
dA). 

The proof that (2.1) obeys (1.1) is as follows (see Fig. 1) [1]. Denote the minimizing 
surface and region m{A) and r[A) respectively. We now consider the region r{ABC) := 
r[AB) U r{BC), and decompose its boundary into the part lying along the boundary of the 
spacetime, which is AB U BC = ABC, and the part lying in the bulk, which we denote 
m{ABC): dr{ABC) = ABC[Jm{ABC), i.e. m{ABC) ~ ABC. Hence 

S{ABC) = area(m(A5C)) < area(m(^SC7)) . (2.2) 

Similarly defining rh{B) as the part of d{r{AB) D r{BC)) lying in the bulk, we have 

S{B) = ^ area(m(S)) < ^ area(m(S)) . (2.3) 

The last step is to note that 

rh{ABC) U m{B) = m{AB) U m{BC) (2.4) 

(this should be clear from the sketch in Fig. 1, but can easily be shown rigorously by decom- 
posing m(AB) into the part lying inside r{BC) and the part lying outside it, and similarly 
with m{BC)), so that 

-— — aTea(m(ABC)) + — — area(m(i?)) = -— — aTea(m(AB)) + ——— aTea(m(BC)) 

4Gn 4Gn 4Gn 4Gn (2.5) 

= S{AB) + S{BC) . 



The area is evaluated with respect to the Einstein-frame metric. No other fields enter into the formula. 
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ABC ABC 




Figure 1: Sketch of the different regions and surfaces appearing in the holographic proof of SSA. 
A, B, C are boundary regions. Each region, as weU as their unions AB, BC, and ABC, has a corre- 
sponding bulk region r{AB) etc. The bulk region f{ABC) is defined as r{AB)Ur{BC) (left diagram), 
while f{B) — r{AB) n r{BC) (right diagram). Equations (2.2), (2.3), and (2.4) are hopefully clear 
from the diagrams. 

Combining (2.2), (2.3), and (2.5), we arrive at (1.1). A similar argument leads to (1.2).^ 

This proof uses precisely the input that is given by the RT formula, namely the homology 
and minimal- area conditions. Nothing more needs to be assumed about the topology or 
geometry of the bulk spacetime. (Not even the null energy condition needs to be imposed.) 
Such generality is both a blessing and a curse, as no constraints on the spacetime can be 
deduced by requiring that SSA be satisfied. In fact, even the minimal-area condition can be 
relaxed a bit, as the proof goes through if we replace the area functional in (2.1) with any 
extensive functional of m [1]. Indeed, it is believed that, in the presence of higher-derivative 
corrections to the classical bulk gravitational action, such as those arising in string theory, 
the area functional gets corrected [5, 18] (although the form of the latter corrections is known 
only in certain cases [19, 20]); if this is true, then SSA is automatically safe, and we learn 
nothing about the form of the corrections (aside from their extensivity) . The same does not 
apply to bulk quantum corrections, whose effect on the entanglement entropy is not simply 
to correct the functional being minimized in (2.1), as can be argued for example by looking 

®It was recently shown [17] using a similar method that the RT formula also implies the inequality 

S{AB) + S{BC) + S{AC) - S{ABC) - S{A) - S{B) - 5(C) > . (2.6) 

Unlike SSA, however, this inequality is not obeyed by general quantum systems, but rather is a special property 
of holographic theories. 
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at the mutual information between distant regions [18]; otherwise httle is known about these 
corrections, whose form may be usefully constrained by SSA. 

2.2 HRT 

It is clear that a major limitation of the RT formula is the restriction to static bulk spacetimes 
and to boundary regions sitting inside constant-time slices. If we search for a covariant 
generalization, then we can no longer appeal to the constant-time slice in the bulk, with its 
Euclidean induced geometry, but rather have to contend with the full Lorentzian spacetime. 
In a Lorentzian spacetime, the area of a spacelike surface with fixed boundary conditions can 
be made arbitrarily close to zero by adding wiggles in the time direction, so the minimal-area 
prescription does not make sense and the RT formula must be altered. Motivated in part by 
Bousso's covariant entropy bound [13, 14], Hubeny, Rangamani, and Takayanagi (HRT) [4] 
suggested replacing the minimal surface with an extremal one, i.e. a spacelike surface that 
extremizes the area functional, or equivalently has vanishing mean curvature; in case there 
are multiple extremal surfaces, we are instructed to choose the one with the smallest area. 
So we have 

^(^) ~ Ar< ™™ area(m) . (2-7) 

4GrN extremal m^A 

(As in the RT formula, m has the same dimensionality as A, so it is codimension two.) The 
choice to replace a minimization procedure by an extremization one when passing from a 
static situation to a dynamical one is a very natural one; for example in classical mechanics 
we go from minimizing the energy for a static solution to extremizing the action for a time- 
dependent one. Similarly, while there is no path-integral derivation of either the RT or HRT 
formula, the HRT formula is reminiscent of the prescription of Euclidean quantum gravity; 
there, since the Euclidean Einstein-Hilbert action is unbounded below, rather than minimizing 
it, one looks for the extremum with the smallest action. 

The HRT formula passes several basic consistency checks. For example, it correctly 
assigns the same entropy to any two regions with the same (boundary) causal domain.^ It 
also correctly reproduces the entanglement entropy in a few simple examples where the latter 
is known, such as for a two-dimensional CFT in at finite temperature [4]. Nonetheless, it is 
fair to say that the HRT formula has so far been subjected to far fewer tests than the RT 
formula. In particular, a crucial consistency test is whether it obeys SSA. Unfortunately, 
the proof of SSA given above does not admit a straightforward generalization to the HRT 
formula. Firstly, because the regions r{AB),r{BC) are not restricted to a given constant- 
time slice, one cannot necessarily take their union and intersection to obtain sensible spacelike 
regions. Secondly, m{B) and m{ABC) are not minimal surfaces, but merely minimal among 

'^In the field theory, this property follows from the fact that the respective reduced density matrices 
p{A), p{A') are related by a unitary transformation. Using the HRT formula, it can be argued as follows. 
The region A is clearly homologous to the future boundary of its causal domain. A' is as well, and since they 
have the same causal domain, A ~ A' . Hence in (2.7) the right-hand side is minimized over the same set of 
surfaces m. 
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the extremal surfaces, so even if one somehow defined appropriate surfaces m{ABC),rh(B), 
one would not automatically have (2.2) and (2.3). 

3. Constant-time intervals 

In this paper we consider a two-dimensional conformal field theory that admits a holographi- 
cally dual description as classical three-dimensional Einstein gravity. Because we will employ 
a Vaidya solution, the bulk theory should admit some form of matter whose stress tensor is 
(at least in some limit) that of null dust; the theory is otherwise unspecified. The field theory 
will be on Minkowski space (with time and space coordinates denoted x,t respectively), and 
we will consider two types of bulk solutions: the planar BTZ black hole, and the planar thin- 
shell AdSs- Vaidya solution (with infalling matter). The corresponding states (or processes) 
in the field theory break both its conformal and Lorentz symmetry, but preserve its spatial 
translation symmetry. The black hole preserves the time-translation symmetry, while the 
Vaidya solution breaks it. 

In general, the entanglement entropy of a spatial region depends only on its causal domain. 
In two-dimensional Minkowski space, the causal domain of a connected region (i.e. a spacelike 
curve) is determined by its endpoints; hence we need not concern ourselves with the shape of 
the curve itself. Since we consider states that are translationally invariant in x, the entropy 
is a function of the times ti,t2 of the endpoints and of their spatial separation Ax. In this 
section we will restrict our attention to intervals that sit at constant time, so ti = t2 = t. In 
subsection 3.1 we review the fact that the entropy of such intervals obeys SSA if and only if it 
is a concave and monotone- increasing function of Ax (see for example [21] for a more detailed 
discussion). This is a very useful criterion, allowing us to see at a glance in any particular 
example whether SSA is satisfied, while also helping to build our intuition about the meaning 
of SSA. In subsection 3.2 we will study the entanglement entropies of constant-time intervals 
in the Vaidya spacetime, and show that they obey this condition; most of the analysis will be 
exact, but since the entropy is obtained as a rather complicated implicit function of Ax, we 
give a proof-by-plot that the function is indeed increasing and concave. Finally, in order to 
investigate the possible role of the null energy condition (NEC) in enforcing SSA, in subsection 
3.3 we consider the case of a Vaidya spacetime with infalling negative-energy matter, which 
makes a transition from a black- hole to an empty AdS spacetime. We find indeed that, while 
the entropy is still an increasing function of Ax in this case, it is no longer concave (so (1.2) 
is satisfied, while (1.1) is violated). We explain how this result could have been expected by 
comparing the entropy curves for AdS and BTZ. 

3.1 SSA, concavity, and monotonicity 

Consider a system with one infinite spatial dimension, and denote the coordinate x. We 
work at a fixed time t, and suppress the dependence of all entropies on it. If the system 
has translational symmetry in x, then the entropy of a single interval A depends only on its 
length, S{A) = s{lx), where s : M"*" — )■ M. In this subsection we will show that the entropies 
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of all such intervals obey (1.1) if and only if s is concave^, and obey (1.2) if and only if it 
is non-decreasing. As mentioned in the introduction, in order for the regions AB, BC, ABC 
appearing in (1.1), (1.2) also to be single intervals, we assume that A,B,C are adjacent and 
in that order. 

Let a, b, c be the lengths of A, B, C respectively (so -S'(^) = s(a), S{AB) = s{a + b), etc.). 
Set y = c/{a + c). We have 

a + b = yb + {l-y){a + b + c), b + c = {1 - y)b + y{a + b + c) . (3.1) 

If s is concave then, since < y < 1, 

s{a + b)>ys{b) + {l-y)s{a + b + c), s{b + c) > {1 - y)s{b) + ys{a + b + c) . (3.2) 

Adding these inequalities, we obtain (1.1). On the other hand, if s is not concave then there 
exist positive numbers a, b such that s{a + 6) < ^{s{2a + b) + s{b)); choosing c = a we obtain 
s{a + b) + s{b + c) < s{a + b + c) + s{b), violating (1.1). 
Similarly, if s in non-decreasing then clearly 

s{a + b)>s{a), s{b + c) > s{c) . (3.3) 

Adding these inequalities we obtain (1.2). Conversely if s is somewhere decreasing then there 
exist a, b such that s{a + b) < s{a); choosing c = a we arrive at a violation of (1.2). 

The simplest examples are provided by conformal field theories at zero and finite tem- 
perature. The entropy of an interval is given at zero temperature by [22] 

s(y = |ln(^), (3.4) 

and at temperature T > by [23] 

s(/^) = ^ In (^-^ sinh(7rT/^.)) (3.5) 

where c is the theory's central charge and e is an ultraviolet cutoff length. Both functions are 
clearly increasing and concave. For theories with holographic duals, the vacuum and thermal 
states are represented by AdSs and the planar BTZ black hole respectively. Since these are 
static and we are considering constant-time intervals, the RT formula may be applied in order 
to obtain (3. 4), (3. 5) [2, 3], and the SSA property follows from the Headrick-Takayanagi proof. 
We now turn to the AdS-Vaidya spacetime, where this is not the case. 

^Throughout this paper, by concave we mean weakly concave, i.e. for < y < 1, s{yli + (1 — y)l2) > 
ys{h) + {l-y)s{l2). 
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3.2 Vaidya spacetime 



The planar AdSa-Vaidya metric is 



(r^ — m{v))dv'^ + 2drdv + r^dx^ . 



(3.6) 



The boundary is at r — )• cxd, where the v coordinate is identified with the t coordinate. The 
spacetime satisfies the NEC as long as m'{v) > 0. We will consider two thin-shell solutions, 
in which the infalling matter is concentrated in a delta- function at v = 0. In this subsection, 
we will take 



where m > 0, representing a shell of (positive-energy) matter in otherwise empty AdS which 
collapses to form a BTZ black hole. In subsection 3.3, we will consider the opposite situation. 



in which a shell of negative-energy matter collapses onto and annihilates a BTZ black hole, 
leaving behind empty AdS. 

In subsubsection 3.2.1, we construct and discuss the Penrose diagram for the positive- 
energy Vaidya spacetime. This sets the stage for the systematic study in 3.2.2 of the global 
structure of the spacelike geodesies that are symmetric under x — t- — x, since those are the 
ones that give the entanglement entropies of constant-time boundary intervals. We follow 
the method of [6] to solve the geodesic equation in the thin-shell geometry. As noted in 
[7], an interesting feature is that some of the geodesies pass behind the horizon; in other 
words, accordng to HRT the entanglement entropies of some boundary intervals depend on 
the metric in the region behind the horizon. Finally, we use the solutions to test SSA. Given 
the theorem proven in subsection 3.1, this amounts to showing that the length of the geodesic 
is a concave, monotone function of the separation of its endpoints on the boundary. This is 
shown in Fig. 4. On the other hand, it fails to be concave in the negative-energy case, as 
shown in Fig. 6. A heuristic way to understand the reason for this behavior is the following: 
A geodesic that is anchored at a time t > will, if the interval is sufficiently short, lie entirely 
in the AdS region. But if we increase the length of the interval, the geodesic will eventually 
cross into the BTZ region. So the entanglement entropy s{lx) will cross over from (3.4) to 
(3.5). As (3.5) has a larger slope than (3.4) for the relevant values of Ix, concavity fails in the 
crossover region. 

3.2.1 Penrose diagram 

Before we solve the geodesies, it is helpful to plot the Penrose diagram of Vaidya space, and 
to have some ideas on the global structure of the spacetime. If the shell is thin, one can 




(3.7) 




(3.8) 
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consider the two regions, pure AdS3 and BTZ separately, and match them together along the 
shell later on. 

If we define /(r, v) := — m{v), and introduce 

dv 

v = t + f, •with df = —. r, (3.9) 

f{r, V) 

one gets the familiar metric, 

' -r^dt"^ + K + r'^dx^, iiv<0, AdSs 



ds^ 



^ ^2 



[-{r^- m)dt^ + + r^dx\ if > 0, BTZ 
One can easily solve, for AdSs, 

f = -l/r (3.10) 

while for BTZ, 



1 a/ Tfi 

arctanhf^^ — ]. (3.11) 
m r 



As we will discuss, the mass parameter m is free to set to 1, however if we want to simulate 
a thick shell spacetime, a variable m would be necessary. Note that for x > 1, 

TT 

arctanh(x) = arctanh(l/a;) — i — . (3-12) 

To find the Penrose diagram of AdSa, the standard method is to introduce u = t — r, and 
apply a series of coordinate transformations, 

j V = t + f j V = ln V j V = tanV 

yu = t-r \^u = - ln{-U) yU = tanU 

For BTZ, one should better introduce some coefficients, for example, 

1 V -Vn 

v= - In V, F = 6 tan -. (3.13) 

a c 

However it turns out only a is necessary, and the real value requires a = ^fm. Thus one has 

{v = t^f {v=^\nV fF = tany 

\u = t-r |u = -^ln(-?7) |i7 = tanC/ 

Since (r, w) are global coordinates, the (V^, f/) should be matched at the shell. With 
2f = V — u, the coordinate matching at v = shell tells us 

UA 



The matching in [/ is, 

2 



tanh ( ) . (3.14) 



ln[tan(-?7A)] = ^ tanh ^ln(-tanC7B) , (3.15) 
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or 



Ub 



arctan 



2 arctanh 



ln(— tan Ua) 



(3.16) 



The Penrose diagram of a thin-shell Vaidya space is given in Fig. 2.^ The asymptotic AdS 
boundary is intentionally kept as vertical, while the behavior of the singularity is more tricky. 
With the nontrivial matching of coordinates along the shell, it seems that we are not able to 
find a horizontal straight singularity in Penrose diagram. 




Figure 2: Penrose diagram for a thin-shell positive energy Vaidya space. The 45-degree red line is the 
infalling null shell. The vertical blue line is the asymptotic AdS boundary, and the black curve on the 
top is r = singularity. The dashed 45-degrce lines are BTZ's horizons and their extension back to AdS 
region. The dotted curves indicate constant radius tracks. A series of coordinate transformations are 
imposed on BTZ region, to match the coordinates in AdS along the mass-shell. The asymptotic AdS 
boundary of BTZ spacetime is intentionally set to be vertical. However the black curved singularity 
r = is in general not able to keep horizontal straight simultaneously. 



To mimic a thick-shell Vaidya space, one could impose multiple mass shells in the space, each of which 
emitted from the asymptotic boundary at different time. The matching of coordinates is much more messy. 
Keeping both singularity and boundary straight simultaneously is a very tricky work, and in general we believe 
it is not possible. Nevertheless it provides us an analytic method for general constructions. 
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3.2.2 Spacelike geodesies 

To solve the spacelike geodesies, the equations of motion could be found, 

- f{r,v)v'^ + 2r'v' + = ^, 



Pi 



2 2/2 

r — r V 



rv" + 2r'v' = 0, (3.17) 



where the conserved momentum px = r x is introduced. In this note, the superscript prime 
indicates and the overdot means the derivative by the affine parameter r. Since the 
metric components do not contain t except for along the thin-shell, it would be helpful to 
introduce another "conserved" momentum along t, 



E = fir,v)i 



r'^i, if f < 

(r^ — m)i, if f > 



The system has a scaling invariance under Vh = ^frn^ 

{r^Vhr, t^—, x^—, v^—, Px-^ThPx, E^rhE}. (3.18) 
rh Th Th 

Equivalently one is free to set the mass m = 1. It makes sense since in Poincare coordinates, 
there is only one dimensionful scale. 

For a spacelike interval on the boundary, time t = constant^ and due to the symmetry 
of the spacetime, the geodesic anchored at the endpoints of the interval is symmetric. This 
simplifies our calculation. 

For simplicity, the thin shell is set along u = 0. On the other hand, since both AdS and 
BTZ has f(r oo) = 0, then the endpoints of geodesic have 

"^b = "^boundary ~ (^ "I" 'f'^boundary — tb- (3.19) 

If th < 0, the whole spacelike geodesic is in AdS bulk, while > 0, it might cross the shell. 
For a spacetime with metric, 

ds'^ = -f{r)dt^ + + r'^dx'^, (3.20) 

since there is no x-dependence and i-dependence in metric, two momenta are conserved, 

Pt = -E = gttP^ = -f{r)i, p^ = gxxp" = r^x. (3.21) 
The affine parameter condition is 

dx" dx" ^ 

which leads us 



f' = f + E^- f,pl, {r'f = "^ + "^,E-- r'f. (3-23) 

Px Px 
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One can solve the minimum of radius r with r(r^,) = r'(r^:) = 0. For a geodesic with constant 
time, i.e. E = 0, one finds r^: = p^. 

If the geodesic is totahy in AdS^ space, the symmetry requires E = 0, and 



•2 2 2 / /\2 6/2 4 

r = r — p^, [r ) = r /p^ — r , = p^. 



(3.24) 



One can calculate the proper length L and the expansion / along x, 



LAdS =JdT = Jdr/r = 2i: 
lAdS =Idx = J dr/r' = 2 /; 



dr 



21og(r + V"^ Pz) = 21og(2roo) - 21ogpa;, 



dr 



r2. /rl_i 
pi 



2 

Pa;' 



It is clear that 



LAdS = 21og(2roo) + 21og 



lAdS 



Similarly in BTZ, with E = 0, and set m = 1, 



^2 _ / 2 



„'2 



l)r2 



PS 



Pa 



(3.25) 



(3.26) 



(3.27) 



Lbtz = 21og(2roo) - 2 log ^^1^, 

Ibtz = 2 arctanh(p2:) + ivr. (3.28) 

One could find the relation, 

Lbtz = 2 log(2roo) + 2 log (^sinh (^^) ) ■ (3.29) 

A crucial feature is Lbtz{1) > LAds{l) for any I > 0. 

If we shift a pure spacelike interval on the boundary, with fixed length /, along time 
direction from AdS region to BTZ region, when tf, < 0, the corresponding geodesic is totally 
in AdS bulk. After it passes over v = 0, the geodesic crosses the shell for twice. And if i;, > |, 
the geodesic will be totally in BTZ bulk. Since the ones totally in AdS or in BTZ bulk are 
pretty simple, we should focus on the ones crossing the shell below. 

In general, a symmetric spacelike geodesic crossing the shell is nontrivial. Its turning point 
is in AdS region, and after the shell the geodesic might probe region behind the horizon, and 
even falls to the singularity. In Fig. 3 several geodesies are plotted. Any symmetric spacelike 
geodesic with its turning point in the shaded region will falls to the singularity. Note after 
the infalling mass shell, the shaded region is defined by the horizon. 

A geodesic with both endpoints anchored to a pure spacelike interval on the boundary 
has two free parameters, {tb = Vb,l)- Suppose it crosses the shell v = at r = rc, together 
with the conserved momentum px = r'^x, the system can be described by the parameter 
combination {px,rc) as well. 
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Figure 3: Several geodesics(green curves) are shown in tiie figure. Any symmetric geodesic with 
the symmetric point in the shaded region falls to the singularity, while all the others will eventually 
approach the asymptotic AdS boundary. The constant time (orange) curve through the crossing point 
of the horizon and mass-shell is drawn for comparison. The second geodesic from top is interesting, 
that it passes through behind horizon region and reaches asymptotic AdS boundary. Some numeric 
detail is vq — —0.4, and Px = {0.5,0.65,0.75,0.95} respectively for the geodesies from top to bottom. 



Since the geodesic is symmetric, the part in AdS region has Ea = r'^i = 0, i.e. it is along 
a constant time curve. It is easy to find the position of the turning point, (ro,t'o), 

111 , , 

ro=Px, vo = to = . (3.30) 

^'o rc Px 

It is not hard to solve the geodesic equations, 

r=le-^{e''+pl). (3.31) 
It crosses the shell at r = rc, with affine parameter value at the two crossing positions 

ri,^) =log(reTxAf^). (3.32) 
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The used as superscript and subscript later indicates AdS, while "i?" for BTZ to clear 
up the confusion. One can calculate the expansion along x and the proper length, 



Tc2 



Tel fcPx 



La = r. 



{A) JA) _ ( + Vrg - pI 
; - \/rl-pl 



c2 'cl 



log 



(3.33) 
(3.34) 



Since the geodesic is symmetric, one can focus on the half part with r increasing from r = rp 
to the boundary. Then at the interception with the shell. 



/ 



rt. 



= —^/'^I-pI, r's 

rc Px 



1 



E^EB = -:^^^Jrl-pl. (3.35) 



E ^ and the geodesic experiences a refraction. 

Now let us turn to the BTZ part. With a nonzero E, the solution is 

i?i = r2 = ^ (e" + B+e-^) (e" + 5_e-") = ^ (« + B+) (a + B. 



where a = e^"^, and 



For simplicity one can introduce 



B± = (p, ± 1)2 - E' 



(3.36) 
(3.37) 



A±=pl-{l±Ef, 



(3.38) 



and they satisfy 



^(i?+ + S_)-l = J(^+ + A_) 



It is easy to verify 



1 

4 



rlrV - 1 = - (e^ + ^+6"^) (e^ + A.e"^) . 



Some details of a general solution is discussed in section 4.1. 
At the crossing point, 



Or 



1 



-{B+ + B^) + Arl± ^y-AB+B^ + {B+ + B^ - ArlY 



(3.39) 
(3.40) 



(3.41) 



where ± is the sign of 



. One can solve 



AxB = - log 



B+ + ac 
B^ + 



Ats = - log 



+ ac\ , All 

— , Lb = — Tc, tb = Vb = At^+arctanh — 

A+ + acJ Tc 

(3.42) 
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Combining all the three parts together, 




(3.43) 
(3.44) 

We have already regularized the proper length by removing the coimnon infinity part 2 log(2rQQ) 
It might be useful to write the boundary time value as 

e-^ = + (3.45) 

To verify SSA of pure spacelike intervals on boundary, we need to keep fixed, and check 
the relation between L^eg and Ix- 

Unfortunately it is hard to derive a clean formula of L(lx), so we turn to numeric solution. 
In Fig. 4, we give a L(lx) function example with tb = 0.5 on the asymptotic AdS boundary. 
The proper length L increases monotonically with Ix- When Ix is small, it behaves like BTZ 
space, and the large Ix behavior is determined by AdS part, which is reasonable, since small 
expansion geodesic is mainly in BTZ bulk while large expansion ones are mainly in AdS space. 



L 




Figure 4: L{lx) for pure spacelike interval at tf, = {0.5,1.1,2} on the boundary, respectively from 
bottom to top. The black curves are L{lx) functions for BTZ and pure AdS, from top to bottom 
respectively. 



To verify SSA inequalities, one can define 



hiA, B, C) = Sab + Sbc -Sa- Sc, (3.46) 
hiA, B, C) = Sab + Sbc - Sabc - Sb- (3.47) 
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If both (/i,/2) are nonnegative, SSA is satisfied. 

The L{lx) function in Fig. 4 tehs us that it is also concave, and thus SSA is satisfied. In 
addition, one can check numerically of I2 function, as in Fig. 5. We fix {Ia, Ib, Ic) = (0.5, 1, 0.5) 
and move the pure spacelike intervals on boundary from AdS to BTZ. The bump behavior of 
I2 indicates that SSA is always satisfied. "'^'^ 



h 

0.30 r 




Figure 5: With {IaJbJc) = (0.5,1,0.5), we move the pure spacelike intervals on boundary from 
AdS to BTZ. I2 is always positive. The green part is for AdS, and the red part is for the case that all 
the relevant geodesies are in BTZ space. 



3.3 Negative-energy Vaidya spacetime 

In physical Vaidya space, a spacelike geodesic could probe from BTZ region into AdS region, 
and experience a flatter function of L{lx)- This induces us the idea that how the function 
L{lx) should behave if the two bulk regions are exchanged. It turns out the exchange is 
connected with the violation of NEC, and SSA is violated. 
For a general Vaidya space, 

ds^ = - fir, v)dt^ + — r + r'^dx'^ = - fir, v)dv'^ + 2drdv + r'^dx^, (3.48) 

f{r,v) 

where 

f{r,v) = r'^ -m{v). (3.49) 

The NEC requires 

Tf^yuf^n" > 0, (3.50) 

^"Thick-shell Vaidya is studied in [7], the monotonicity and concavity of L{1) are apparent in Fig. 3 of it, 
as well as SSA. 
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where n'^ is lightlike, i.e. n^^n^ = 0. In a Vaidya space with asymptotically AdS^+i boundary, 
the Einstein equation is 

G^,^ + Ag^,^ = T^^, where A = (3-51) 

here I is the AdS radius, and A = —3 for d = 2 if setting / = 1. Suppose a null vector 
is = (a, 6,0), the x component is set to zero for convenience. Solving the null normal 
equation n^n^ = gives 

a = 0, or b = f{r,v)a/2. (3.52) 

The first one is simply just ~ (0, 1,0) exactly along coordinate r, while the second ~ 
(1, //2, 0). The first one is trivial since Grr = T^r = 0, while the second has, 

T77' 

T^.n^n'' = G^.n^'n^ = — . (3.53) 

Thus the NEC is equivalent to m'{v) > 0. 

If we turn on a negative-energy mass shell, with m'{v) < 0, it violates the NEC and 
results into a negative-energy Vaidya spacetime: a BTZ evolves into a pure AdS spacetime. 

For an observer sitting on the asymptotic boundary, if interested only in the pure spacelike 
intervals with expansion along x of Ix, there are basically three kinds of spacelike geodesies. 
As usual, suppose the thin shell is emitted from the boundary at v = 0, and the observer is 
at t = tb- 

1. If tb < 0, the geodesies is totally in BTZ region. 

2. For tb > 0, if Ix < lx{tb), it is totally in AdS^. 

3. For tb > 0, if Ix > Ixitb), the geodesic crosses the mass shell, and approaches BTZ 
region. 

It is easy to find lx{tb) = 2tb, which is the same as in physical Vaidya space. 

We shall keep the details in appendix A, and list some features here. The calculation is 
complicated, and the geodesic takes different branches of the solution when the pure spacelike 
interval on boundary changes with parameters {tb, Ix)- In Fig. 6, we plot a series of geodesies 
with different boundary time value of t;,. A general behavior is when tb increases, the curve 
moves from BTZ to AdS. From the value ranges of tb, (A. 24, A. 30, A. 37), there are some 
overlaps. The combinations give continuous L{lx) curves. An important property is although 
the curves are still monotonic, the concavity is not preserved, and this results into the violation 
of strong subadditivity. 

Keeping tb > fixed, the spacelike geodesies with small Ix feels only the effect of AdS^ 
bulk, while the ones with large Ix cross the shell, and eventually the BTZ's effect dominates 
if Ix is large enough. 

In both Fig. 6 and Fig. 4, for a pure spacelike interval on the boundary, the small Ix 
behavior is determined by the bulk after the mass shell, while the large Ix is dominated by 
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4 



^ Ix 
12 



Figure 6: A series of geodesies ending on AdS boundary at ti, — {0.5, 1, 1.3, 2, 2.5, 3, 5} from top to 
bottom respectively. Red curves are calculated from Tc > 1, the blue ones are from ^ < < 1, and 
the brown come from < ^. Clearly the middle three with ti, = {1.3,2,2.5} are combination from 
two cases. The top and bottom black curves are given for BTZ and pure AdS^ respectively. 

the one before the shell. Since LsTziP) > LAdsi^)-, only the L{1) curve for the negative-energy 
Vaidya contains a non-concave part, and thus SSA is violated. Heuristically we could make 
a conjecture. Suppose there are two adjacent regions A and B, sharing the same asymptotic 
boundary. For spacelike geodesic with its endpoints attached on the asymptotic boundary, 
Lsil) > La{1)- Then for any such spacelike geodesic probing from bulk region A into B, SSA 
is violated. On the other hand, if the geodesic is probing from B into A, SSA is satisfied. 
A manifest possible example is a domain wall at constant radius, the two sides of which 
might have very different behavior, say one is AdS while the other is BTZ, or one is AdS and 
the other is Lifshitz-like space, or the two have different cosmological constants. Naively to 
think, the violation of SSA is determined by properties of the two regions, while the domain 
wall might impose some effect, e.g. a negative tension might change the story from with a 
positive tension. (This may be related to the fact that proposed entanglement-entropy-based 
holographic c-theorems make use of the NEC, albeit in a static context [15, 16].) However, 
to make the system stable, NEC might have to be satisfied. 

4. General spacelike intervals 

In this section we expand our considerations to include spacelike intervals that do not lie in a 
constant-time surface (with respect to the Lorentz frame distinguished by the system's spatial 
translation symmetry), a case that has been somewhat neglected in the literature so far. In 
particular, in order to check SSA thoroughly we consider a variety of possible configurations, 
in which the intervals are not only not constant-time but not collinear with each other (see 
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Fig. 10). 

It is instructive first to consider the case of general spacelike intervals in a CFT at finite 
temperature. Since the state is translationally invariant in both space and time, the entropy 
is a function of Ax and At.^^ Defining null coordinates u = x — t, v = x + t, the entropy can 
be expressed in terms of the function s defined in (3.5) [4] 

S = l{s{Au) + s{Av)). (4.1) 

(The same formula applies at zero temperature, with s given by (3.4).) Equation (4.1) clearly 
obeys SSA, since each term s{Au) and s{Av) obeys it separately. It was shown in [4] that 
the HRT formula, applied to the BTZ spacetime, reproduces (4.1). 

Note that equation (3.5) has the following interesting feature (for any concave function 
s). If we consider the quantity I2 defined in (3.46), which measures how far (1.1) is from 
being saturated, it can be written 

i2 = l{i^ + r2), (4.2) 

where 

I2 = s{AuA + Aub) + s{AuB + Auc) - s{Aua + Aub + Auc) - s{Aub) (4.3) 

and similarly for /I- Now, in the limit that A is null, either Au^ — )• or Ava (depending 
on which way it is tilted), implying either I2 ^ 01 I2 ^ 0. If we also take the limit where 
C is null, but tilted the other way — as in the trapezoid configuration of Fig. 10 — then both 
I2 and I| will go to 0, so /2 — )■ and SSA will be saturated. (Note that, even when both A 
and C are null, all of the intervals appearing in the definition of I2 are spacelike.) In the rest 
of this section we will consider a Vaidya spacetime that makes a transition from AdS to BTZ 
(or vice versa for negative-energy Vaidya). We will pay particular attention to the trapezoid 
configuration with A and C close to null: since SSA is just barely satisfied in both AdS and 
BTZ for this configuration, it offers the best chance of finding a violation if one exists. 

As shown in Fig. 7, for a spacelike geodesic with both endpoints attached on the bound- 
ary, there are four possible configurations: it could be totally in AdS or BTZ, or it might 
cross the thin-shell once or twice. We will discuss these possibilities in subsections 4.1 and 
4.2, finding exact solutions in each case. Then in subsection 4.3 we will put together the 
solutions in order to test SSA. This involves a certain amount of numerical work, since the 
formulae for the proper lengths of the geodesies in terms of boundary data are complicated 
and implicit. We find that the infalling matter actually has the effect of increasing I2 for the 
trapezoidal configuration, as long as it has positive mass (see Fig. 13), so SSA is safe; nor 
do we find violations for Ii, or for either Ii or I2 for any other configurations. On the other 
hand, we do find that negative-energy infalling matter does drive I2 negative (see Fig. 15), 
showing that violation of SSA is correlated with violation of the null energy condition (NEC). 

^'^From now on we will use Aa; instead of as in pure spacelike case in previous section, to avoid the 
potential confusion to the length I = ^ (Aa;)^ + (Af)^ on the boundary. 
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shell 




Boundary 



Figure 7: The four basic kinds of geodesic with both endpoints attached on the asymptotic AdS 
boundary. The geodesic might be totally in AdS bulk or BTZ space, or it might crosses the thin-shell 
once or twice. Note that the shell-crossing geodesies may reach the region behind the horizon (not 
shown in the figure). 



4.1 Geodesies totally in pure AdS or BTZ 

The formula in pure AdS or BTZ space is simple. One can easily find, in pure AdS, the 
solution is 



1 



where the conserved momenta are defined 

E = r'^i, Px = r'^x. 
The boundary has the affine parameter values 

Too = In 2roo , T_oo = - In 2roo + ln(p^. - E^) 

while the turning point is at 

2 2 27^2 2 

' turning = Px- ^ <Px 

instead of r = as with E = case, and 

^turning — ^^{Px ^ )• 



(4.4) 
(4.5) 

(4.6) 
(4.7) 

(4.8) 



One could find the time values with {px,E,vo), where vq is the v value of the turning point. 



to = vo + 



1 ''max 



E 



Px ^0 



) '"mm 



to 



E 



Px ^0 



(4.9) 
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and the expansion on boundary, 



^- = ^1^' At = ^, (4.10) 

L = 2\n2roc-ln{pl- E^) ~ ln(Ax + At) + ln(Ax - At). (4.11) 



The discussion in BTZ space is shghtly comphcated. For the geodesic equations 

r^ = f + E^-^pl r'^ = ':^ + !:^E^-r^f, E = ft, = r^x, (4.12) 

Pi Pi 

there are two branches of solution, 

Ri{r) = r\T) = \{e^ + B+e-^) (e^ + B^e'^) , (4.13) 

i?2(A) = r{Xf = [e^B+ - e~^) (e^i?. - e^^) , (4.14) 

here we used the function B± introduced in equation (3.37), and hst them here to be clear, 

B± = {p,±lf-E^ A±=pl-{l±Ef. (4.15) 

The solution Ri{t) was derived in [6]. If we introduce a = e^'^ or a = e^'*', the solutions are 
simplified into 

Ri(a) = ^(B+ + a)(B^ + a), R2(a) = -^(aB+ - l)(aB^ - 1). (4.16) 
Aa 4a 

With the quantities A±, one has 

Ri{t)-1 = ^ {f + ^+e-") {f + A.e~^) , (4.17) 

i?2(A) - 1 = (e^^+ - e"^) (e^^- - e"^) . (4.18) 

Thus the signs of B± tell us whether the geodesic falls to the singularity r = 0, while the 
signs of A± indicate whether it crosses the horizon r = 1. 

An analysis of the solutions R\^2 by quantities A±, B± is illustrative. The solution Ri{t) 
has at least one limit r — )■ ±oo attached to asymptotic AdS boundary, and can be classified 
into four different kinds(Fig.8): 

1. If A±, B± > 0, the geodesic is totally outside of the horizon(Fig.8(a)). 

2. If A^A^ = Bj^B^ < 0, the geodesic crosses the horizon and reaches the singularity(Fig.8(b)). 

3. If A± < 0, B± > 0, the geodesic crosses the horizon twice, and goes back to AdS 
boundary(Fig.8(c)). 

4. IfA-t, B± < 0, the geodesic crosses the horizon twice and hits the singularity twice(8(d)). 
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(a) = 1.5,S= 1 (b) = 1.2,£; = 0.1 (c) = 0.1, £ = 0.5 (d) = 0.5, £; = 2 

Figure 8: The four cases of the solution Ri{a). The dashed hne is the horizon. 

On the other hand, the solution i?2(«) describes the geodesies falling to the singularity(Fig.9): 

1. = B^B^ < 0, the geodesic crosses the horizon once, and the two endpoints 
approach the AdS boundary and the singularity respectively(Fig.9(a)). 

2. B± > 0, A± < 0, the geodesic is totally inside the horizon and both endpoints fall to 
the singularity(Fig.9(b)). 

3. B± > 0, A± > 0, part of the geodesic is outside of the horizon, while both endpoints 
fall to the singularity(Fig.9(c)). 

4. B± < 0, A± < 0, it gives that the value of R2 is always negative, which is non-physical, 
and not shown in Fig. 9. 






(a) p, = 1.1, £ = 0.2 (b) p, = 0.7, £ = 0.2 (c) p, = 1.9, £ = 0.8 

Figure 9: The three physical cases of the solution R2{a). The dashed line is the horizon. 

For a spacelike geodesic totally in BTZ space with both endpoints on boundary, it can 
take only the branch of solution Ri, with B± > 0, which is obvious from Fig. 8. 
Again it might be useful to write down the affine parameter values. 



Tturning = \ (B+B^) , Too = In (2roo) , T_ 



ln(2roo) + ^ln(S+S_^ 



In addition, the turning point is at 



ro 



to = vq + arctanh 



(4.19) 



(4.20) 
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The interval on boundary has 



tmax 



tmin 

And the expansion is 

Ax = ^ln|±, At=iln^, (4.23) 

L = 21n(2roo)-^ln(S+S_) ~ In (sinh (^3^)) + In (sinh (^^)) (4.24) 
4.2 Geodesies crossing the shell 

For a general spacelike geodesic crossing the shell, the solution has to satisfy the equations 
of motion (3.17). Since the geodesic might reach the region behind BTZ's horizon, it is hard 
to tell which branch of solution, Ri or R2 it takes at a given radius r. One should be careful 
to check the evolution of geodesic with the affine parameter. In the following, we argue that 
for the geodesies we are interested, branch Ri is complete to describe them. With the affine 
parameter transformation, the only section of R2 connecting to asymptotic boundary is 1-1 
mapping to the section of -Ri with B^B^ < 0. In addition, the discussion on 2nd derivative 
below justifies the Ri solution. 
Define a = e^"^, one finds that 

Ri{a) = ^ {a + B+) {a + B^) ^ i R'lia) 

R2{a) = -^{aB+-l){aB.-l) [^{0) 

For solution only the section with A^A- = B^B^ < gives negative i?'/, while all 
sections of R2 have R2 < 0. Since we are only interested in the geodesies with both endpoints 
on AdS boundary, the part in AdS must have r"(r) > 0, i.e. if introducing {3 = , 

r^^^ = Ypi^'+pl-El), =^ % = '^^>^ =^ P">^A. (4.25) 
On the thin-shell, 

M,. ^ rifr2>l,r- >0 ^ i?i 
^ V 2r2y ^ [ If r2 < l,r^ < 0,and > 0,r^ < 0. 

Thus we say Ri is enough to describe the part of geodesic in BTZ region. 
One can derive the following matching equations. 



to 



to 



2 V^- + /B^ 



1^^/ A+ + V5T^ 
2 ^\A- + ^B^ 



2 V^- 



(4.21) 
(4.22) 



_ 1 B+-B- 



v' = 2(r^ - t'b), 



(4.26) 



and in AdS, 
while in BTZ, 
and this gives 
i.e. Eb < Ea- 



y' = t' + -4 = ^ + 4 =^r'B= (l-^]r'A-^, (4.27) 



_ 1 _ 1 — 1 

= - If - (4.29) 



4.2.1 Crossing the shell once 

For a geodesic with both endpoints attached to AdS boundary, and crossing the thin-shell 
once, the general behavior is that it starts from its endpoint, falls into AdS bulk, turns back 
at its turning point, then crosses the thin-shell. After the thin-shcll, if r'^lr^ > 0, the geodesic 
goes straightly to the AdS boundary, otherwise it falls back to the singularity but turns back 
at some place, then goes to the AdS boundary. The sign of r'^\rc is determined by {px, Ea, Tc). 
We list below some details of the simplest case, with r'^\rc > 0. The part in AdS has, 

ri^^ = In (r, + ^rl + El-pl^ , ri^i = In {r^ - ^rl, + E\-pi^ , (4.30) 

and 

= — r^^ , + (4.31) 

t'^' = + -4^, (4.32) 

L(A) = ,iA) _ ,04) . (4.33) 



The part in BTZ has 



and 



-{B+ + B_) + Arl ± ^-4B+B_ + {B+ + B_ - Arlf , (4.34) 
T^^)=ln(2roo). (4.35) 



Some of the other quantities are 

The quantities A±,B± are defined by {px,Eb)- 
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The total expansion on boundary has 
and the total proper length 



(4.37) 
(4.38) 



4.2.2 Crossing the shell twice 

Suppose the geodesic crosses the shell at r = rd and r = rc2 for the bottom and top half 
respectively, it is easy to check that rd > rc2, and one should calculate the three parts 
separately. Again it has three free parameters, say {px-,Ea-,vo). 

For the central AdS part, the turning point is at {tq^vq), where 

1 



pI-E% to = vo + 



(4.39) 



The crossing radius {rci,rc2) satisfy 



1 

rci 



to + EA 



1 



ra (tci - ^Jrl-^ + E\-pi) ^ 



= to + EA 



rc2 



The change of t and x are 



pI-e 



A rc2 (rc2 + \lrl2 + El-pi) ^ 



(4.40) 



(4.41) 



rc2 Tcl 



Ax A = At A * 



Px 

Ea 



and the proper length 
La = In 



rc2 + \lrl^ + E\-pl 



In 



At the crossing r = rd, r < 0, thus 



rA 



'cl 
Px 



rli + E\-pl, 



Ebi = Ea[1 



rli+El-pl 



Px 



2r2 



fA- 



(4.42) 



(4.43) 



(4.44) 



To write the solution of part BTZl, one needs to check the values of B^^^ = B±{px, Ebi) 



IfB^^^B^y < 0, 



r2(r) = \{e-^ + B^^'^e^){e-^ + b'l^c^), 



(4.45) 



with r_oo = - ln(2roo). Too = ln(2roo) - | In (^B^+^B^^^ 
If B^^^B^y > 0, r2(r) = + B^^^e-^){e'- + B^e-^), (4.46) 

with r_oo = - ln(2roo) + ^ In (^J^^i^^) , Too = ln(2roo). 
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More specifically, if B^_}^ B^^ < 0, one can introduce 



2rS 



B^2^ + si^) 



+ 4rl, - ^-AB^^^B^y + [B^_i^+B^y-4rl,y 



and 



1 + 

1 + Si'^aci 



Atfii = - In 



1 + AL^^OcI 

1 + ^i^^aci 



LbI = T^l- 



(4.47) 
(4.48) 



While if B^l^B^}^ > 0, one can introduce 



Oici = e 
and 



2tS 



ijW + ^(1) j + 4^2^ _ J-4B^_^)B^y + + B^y - 4r; 



At 



Bl 



In 



(1) 



^ B 



(1)^ 
(1) 



Ln. = ri-lln{Bi^By) 



The boundary value of v for the starting point of geodesic is 



vbi = ^Bi + = Aifii — arctanh 



1 

rci 



The top BTZ part is similar. At the crossing point r = rc2, ta > 0, and 



ic2 
Px 



Eb2 = Ea 1 



1 



2r2 



2r 



(4.49) 

(4.50) 

(4.51) 
(4.52) 

(4.53) 
(4.54) 



c2 



Note that in general, Eb2 7^ E^i- The value r'^ 
just after this crossing point. One finds 

1 



c2 



is important, and determines the behavior 



t'b 



c2 

(2) 



1 - 



K2 



Ea 



(4.55) 



< and at least one of B 



(2) 



c2 



One has to check the value of B)^ = B±{px, Eb2) first. If 

is negative, the geodesic will fall to the singularity, otherwise it will eventually goes to the 
boundary. Since we are only interested in the ones approaching the boundary, the solution is 
simplified. Similarly as above, one introduces a quantity 



2tB 

ac2 = e = 



- (Sf + B^) + 44 ± y^-45f5(^) + (sf + 5i^)-4r,2,)^ 



, (4.56) 
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where "it" is the sign of r'g . The other quantities are 



AXB2 = 2 



Mb2 = 2 In 



The boundary v value is 



vb2 = AtB2 + arctanh 



rc2 



Combining all the three parts together, 



Ax 
At 



L 



reg 



VB2 - VBl, 

La + Lbi + Lb2- 



4.3 Tests of SSA 



Lbi 



(4.57) 



(4.58) 



(4.59) 
(4.60) 
(4.61) 



We discuss in previous section about the connection between SSA and the NEC for pure 
spacelike intervals. In this subsection, the discussion is extended to general spacelike intervals, 
and we argue the causality somehow preserves the connection. Some examples are checked 
in Vaidya space, verifying SSA. In the end of the subsection, an example in negative-energy 
Vaidya space shows that SSA is also violated when the NEC is violated. 

For general spacelike intervals, there are basically three kinds of combinations, i.e. trape- 
zoid, zigzag and linear combinations, as in Fig. 10. 



B 




Figure 10: The three kinds of boundary intervals to test strong subadditivity: trapezoid, zigzag and 
linear combinations. 



In all of the three cases, the straight line one is the basic, and we will show some features 
of this first. As discussed, for a straight line interval, SSA is equivalent to the combination 
of concavity and monotonically increasing function of L{1), thus we basically need to show 
the properties of the function. In Fig. 11, some comparisons are given. The Vaidya case is 
much closer to pure AdS than to BTZ, which is reasonable, since although the center of the 
straight line interval is at t = on boundary, the interception of the geodesic to the thin- 
shell has t > 0, and the geodesic is mainly in AdS bulk. Given the same parameters of the 
intervals, the corresponding geodesic in BTZ has larger proper length than in AdS, while the 
one crossing the mass shell is between them. Varying the slope At/ Ax of the interval, one 
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Figure 11: L^eg as function of length I for straiglit line intervals on AdS boundary. In the left figure, 
the bulk is pure AdS space, with different slopes, ^ = {0, 0.5, 0.9} from top to bottom respectively. In 
the right figure, a comparison is given for pure AdS, BTZ and thin-shell Vaidya space. For convenience, 
we set = 0.5 and fix the center of line on the shell, v ~ on boundary. The curves are for BTZ, 
Vaidya, and pure AdS from top to bottom respectively. All the curves here are obviously concave and 
monotonically increasing, and thus satisfy strong subadditivity inequalities. 

finds larger slope gives less proper length. In each given example in the figure, the function 
L{1) is monotonically increasing and concave, therefore SSA is preserved. 

In the following, we test SSA for linear, trapezoid and zigzag combinations of intervals 
respectively. For simplicity, the lengths are fixed, Ia = Ic = 0.5, Ib = 1- The intervals are 
shifted from AdS to BTZ, thus boundary time value is one of natural parameters. The other 
variable is the slope of the interval A, from which we determine the slopes of B and C. 

In Fig. 12, SSA is verified for linear combination of intervals. We vary the slope on (x, t)- 
plane, and find the functions of Ii^2 are similar to pure spacelike interval case. The x-axis 
indicates the time of the center of interval i?.(Note that on the boundary, v = t.) 

For the trapezoid combination, for simplicity, we suppose A, C have the same(with op- 
posite signs)slopes on (t, x)-plane, while B is pure spacelike. Several examples of different 
value of At/Ax are given in Fig. 13. As we previously argued, in either pure AdS or BTZ 
space /2 = when j4, C are lightlike and with opposite signs of At/ Ax. The vanishing of I2 
is approximately confirmed by the red curve which has At/ Ax = 0.98 and intervals A, C are 
near lightlike. Nevertheless I2 does not vanish around the transition from AdS to BTZ. 

Now let us check for the zigzag combination. For simplicity, the intervals are intentionally 
to set that A parallel to C, and AB^ BC are pure spacelike. Several examples are illustrated 
in Fig. 14, and for all of these both Ji^2 > 0, verifying SSA. We also find numerical evidence 
that /i,/2 generally do not vanish even if all the three intervals in zigzag combination are 
null. 

Although in both Fig. 13 and Fig. 14, the Ii functions look pretty flat, a close look shows 
that they behave more like a smoothed step function, similar to the left figure of Fig. 12. 

For all the examples examined above, Ii is in general a smoothed step up function, always 
positive, while I2 has a bump in the transition from AdS to BTZ, and there is no sign to 
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-0.4 -0.2 0.2 0.4 0.6 0.8 1.0 -0.5 f 0.5 1. 

Figure 12: SSA is verified witli Ii^2 functions for linear combination of intervals. In the left figure, 
At/ Ax — {0, 0.204, 0.314, 0.436, 0.577} respectively for the curves from bottom to top. In both figures, 
the same color corresponds to the same value of At/ Ax. The length of intervals are set as Ia — Ic = 
0.5,1b = 1. The sharp zigzags are due to numeric errors. The x-axis is the boundary time value of 
interval B. 




Figure 13: SSA is verified with /i_2 functions for trapezoid combination of intervals. In the left 
figure. At/ Ax — {0,0.314,0.75,0.98} respectively for the curves from bottom to top. In both figures, 
the same color corresponds to the same value of At /Ax. The x-axis is the boundary time value of 
interval B. 



negative value region. Therefore one can conclude that both Ii^2 are kept positive and SSA 
is preserved. 

On the other hand, as conjectured, the negative-energy Vaidya space, which violates 
the NEC, should violate the SSA. An explicit example is shown in Fig. 15. Similar to pure 
spacelike interval on the boundary of negative-energy Vaidya bulk, which has a monotonically 
increasing but non-concave L{lx) function(Fig.6), and thus /i > while I2 is not always 
positive. Compared with Fig. 13, for the same combination of intervals on boundary, the Ii^2 
functions seems to exchange the future with the past, and flip the shape upside down, e.g. 
the bumped I2 in Vaidya v.s. the well-shaped I2 function in negative-energy Vaidya space. 
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Figure 14: SSA is verified witii /i^2 functions for zigzag combination of intervals. In the left figure, 
At/ Ax — {0,0.577,0.75,0.98} respectively for the curves from bottom to top. In both figures, the 
same color corresponds to the same value of At /Ax. The x-axis is the boundary time value of the 
centers of the intervals A, B and C. 




Figure 15: The /1.2 functions for negative-energy Vaidya, with slope At/ Ax = 0.7 for interval A. As 
usual, Ia — Ic — 0.5, Ib = 1- The x-axis is the boundary time value of interval B. 



Without time consuming numerical calculation, one can safely say, if the lengths are kept 
Ia = Ic = 0-5, Ib = 1, the trapezoid combination is the easiest to violate SSA, especially 
when A, C are near-null. However, the appearance of negative I2 might depend on the length 
combination one choose. 

In addition, for given sets of parameters, /i_2 are constant in both pure AdS or BTZ 
bulk, so as long as these functions start to change, some of the relevant geodesies crosses the 
thin shell. Since I2 in Fig. 15 does not jump to negative value but varies smoothly, it suggests 
us that violation of NEC, equivalently with m'(v) < in this example, is not a sufficient 
condition for SSA's violation. We would like to prove that more formally in our future work. 
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A. Negative-energy Vaidya: Details of calculation 

Since the constant time geodesies in pure AdS^ or BTZ is discussed already, one needs to 
focus on the ones crossing the thin shell. Starting from the symmetric point (r, v) = {px, vq), 
the geodesic should have constant time behavior in BTZ region, while tAdS after the mass 
shell. Thus it crosses the mass shell at {r,v) = (rc,0). In BTZ, as given by the standard 
form, 



here we have already scaled the mass of the shell m to 1. The constant time behavior in BTZ 
indicates 



From the equations of motion (3.17), one can derive some of the quantities relevant. To avoid 
the Delta function behavior of v" , v' has to be continuous, thus 



here the subscript — and + indicate before and after the mass shell, i.e. BTZ and AdS^, 
respectively. In BTZ, Ebtz = 0, 




(A.l) 




(A.2) 



v'{v = 0) = 2(r+ - r'_) 



(A.3) 




r. 




1) 



v'{v = 0) 




(A.4) 



These together gives 




(A.5) 



and r'j^ reverses its sign aX = On the other hand, in AdS 




(A.6) 



and thus 




(A.7) 



from which one calculates the conserved momentum along t 




(A.8) 
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For the geodesies connecting to AdS boundary, the proper length and expansion along 
x-direction are calculated, 



L 



BTZ 



2 In 



1 + 



I BTZ = 2ArcTanh p. 




In 



Px 



pI 



PxVr 



Px 



LAdS = 21n (r + ^r2-p2+^2 



lAdS 



{pI — E'^)r 



(A.9) 



(A.IO) 



(A.ll) 



(A.12) 



As we discussed, in BTZ there are only two cases for connected symmetric spacelike 
geodesies with E = one is totally outside the horizon, and the other totally behind the 
horizon. In the first case, the radius of geodesic after crossing the shell continuously grows, and 
approaches the AdS boundary eventually. The geodesies behind the horizon has complicated 
behavior, and we will discuss more details below. A close study shows there are four situations, 
as in Fig. 16: 




Figure 16: The four cases of geodesic in negative-energy thin-shell Vaidya space. The vertical straight 
line indicates the thin shell at u = 0, and the horizontal dashed line at r = 1 is the horizon. From top 
to bottom, > 1, ^ < < 1, < i with — E'^ > 0, and < | with pi — E'^ < respectively. 
The geodesies generally encounters a refraction when crossing the shell. The figure is not drawn in 
scale, only showing the qualitative behavior. 



1. Tc > 1, the part in BTZ region is totally outside the horizon. 
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Obviously, it has some features 

rc >^?x > > ^ r+ > 

and 



With the equation 



one can solve 



_|_ _g'2 _ 



r=^e-^{e^^+pl-E^), or r = ln(r + V^^T^E^^^). 
At r = Too— >-oo,r = Too — >■ oo. The geodesies crosses the shell at 



Tc = ln(rc + ^rl + - pi) 

Thus the proper length 

L = J dT = Too-Tc = ln(2roo) - ln(r + ^r^ + E'^ - pi) 

Similarly, 



Px 



rc{rc + ^Jrl + E'^ - pi) 
The total proper length and expansion along x are 

L = Lbtz + '^L, Ix = Ibtz + '^{x — Xc). 



(A.13) 
(A.14) 
(A.15) 

(A.16) 
(A.17) 

(A.18) 
(A.19) 

(A.20) 



One simple way to check the transitions is keeping vb fixed, and change the expansion l^^ 
furthermore, one is able to discuss subadditivity and strong subadditivity properties. 
However in our method, the initial parameters are {Px,vq) of the symmetric point, 
which could be easily transferred into {px^fc)-, and it is in general messy to find (tb,lx) 
numerically. It is helpful if we can calculate the constraint with = thoundary fixed. In 
AdS, 

v' = t' + 



V = t + f = t ■ 



1 , , r' i r' E r' 



+ 

x 



+ 



Px 



(A.21) 



and 



/•boundary 

Vb = + / —dr 

J shell 

oo 

1 E 

H X 

f Px 

shell 



boundary ^/ 



dr 

oo 



rboun 
J shell 



dr 



boundary / IE 



r ( 

shell \ 



2 + / 

r r Pa 



dr 



E 



rc rc{rc + y/rj + E"^ - pi) ' 



(A.22) 
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With tf, fixed, one can solve 



{2rc + tb- 2rlhf 
The constraints from r^. > Px > ^ and tfe > gives that 



^,^ rg(4-4rA + «;) ^^^^^^ 



1 4r 4 

-<t6<^^,, and tfeG(0,-). (A.24) 



For reference, we hst some details of the constraints here. 



„2 



> ^ tfe < 2(r - v^r2 - 1) or > 2(r + v^r2 - 1), or < 



4 + t' 



2 



Px — - - — — V v c -/ ~- — \- ' v c -/7-- L ^, 

4r 

p. > 1 ^ h< (A.25) 

All the above constraints combined together gives us the constraint (A.24). One can 
further write down the formula of {Ix, L), 



2(1 - r^tb) - J 4 - Arch + t 



+2- 



i7r + 2V4-4rci;, + t^ + ln 

2(l-ret6) + -^4-4retb + tg 

L = In f^^(2^c + t.-2.gt,)2N 

V 4r, + tfe - 4r2tf, ^ ^ ^ 

Unfortunately we are not able to find a clean expression of L{lx), which would be very 
helpful for the entanglement entropy discussion. 

Several geodesies ending on AdS boundary at different value of tb are given in Fig. 17. 
Since they are not very far from u = at boundary, the curves are pretty close to 
BTZ curve. For < tf, < 1, the L{lx) curve connects to the bottom AdS curve, which 
happens when the symmetric point in BTZ moves to (f = 0,r = 1). However when 
tb > 1, there is a disconnection, which comes from E > 0. 

2. \ < rl < I, the part in BTZ region is totally behind the horizon, as the second top 
curve in Fig. 16, with some features: 

rc<Px< 1, < 0, r\>Q, E > 0. (A.27) 

Since r(x) is a monotonic increasing function in AdS region, it has 



r = +Vr2 + ^2 _p2^ (A.28) 

and the solution is 

r = ^e-^ (e^^ + pi - E^) , or r = ln(r + y^r^ + E^ - pi). (A.29) 



- 35 ~ 



L 



15 - 



10 - 



5 - 




Figure 17: Three geodesies ending on AdS boundary at tt = (0.5, 1, 1.3) from top to bottom respec- 
tively. Tlie top black curve is for BTZ black hole, and the bottom black one is for pure AdS^. The 
geodesies are pretty close to BTZ curve. 



Most functions are the same as the case Vc > 1, e.g. Tc,L,x,tb,Px- Although ^ < r'^ < 
< 1 here, different from the first case, the equations leads to the same constraint, 

1 4t* 

- <tb< ^^"^ ^ (^-^^^ 

3. < < 2 and pi - E'^ > 0. In this case, 

r'_ < 0, r'^ < 0,£; > 0, (A.31) 

and the sign of r is changed, 

r = -^/r^ + -pI. (A.32) 

The solution takes a slightly different form, 

r = ^e^ [e-^^ + pI - E^) , or r = - ln(r + + E^ - pi). (A.33) 

lipl — E'^ > 0, the radius r falls down first, and turns back at = pl — E'^, approaching 
AdS boundary eventually, at Too = ln(2roo) — ln(p^ — E'^). The geodesic crosses the thin 
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shell at Tc = — ln(rc + y^r^ + E'2 — p'^y One can calculate the expansion and proper 
length 

^ = + vM- m ^ -pI + E') , (A.34) 
L = JdT = Too-Te = ln(2roo) - ln(p^. - E^) + In (r, + ^f^^fT^) . 



(A.35) 



The constraints on are, 



I 3 _ 4j-2 

"^^""^<5^ < < 2r.(l - 2;g) 

when rc = -k, — < t;, < 4, 



r. 



c 



V2 

and 



1 2 

when 5 < rc < ^, — < tf, < — , (A.36) 



t6G(\/2,oo). (A.37) 
4. < < 2 and p^. — E'^K 0. The geodesic falls to the singularity, in finite proper time 



f = -Un{E^-pl). (A.38) 



The proper length is 



/ 



dT = f-Tc = ln(^2 - pI) + In (tc + ^/rl + E^ - pl^ . (A.39) 



Although it reaches the singularity in finite proper time, the expansion along x grows to 
infinity, as well as value of So we argue it is not necessary to consider any reflection 
from the singularity r = and bouncing back to AdS boundary. 
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